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Abstract. We report on the XMM-Newton observation of the HMXRB X Persei, the prototype of the persistent and low- 
luminosity Be/neutron star pulsars, which was performed on February 2003. The source was detected at a luminosity level 
of ~ 1.4 x 10 35 erg s _1 , which is the highest level of the latest three decades. The pulsation period has increased up to 
839.3 s, thus confirming the overall spin-down of the NS detected in the previous observations. The folded light-curve has a 
complex structure, with features not observed at lower luminosities, and shows a significant energy dependence. The spectral 
analysis reveals the presence of a significant excess at low energies over the main power-law spectral component, which can 
be described by a black-body spectrum of high temperature (Web ~ 1.5 keV) and small emitting region (Rbb ~ 340 m); its 
properties are consistent with a polar-cap origin. Phase-resolved spectroscopy shows that the emission spectrum varies along 
the pulse period, but it is not possible to prove whether the thermal component is pulsed or not. 
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1. Introduction 

The X-ray source 4U 0352+309 is the brightest and per- 
haps most famous member of the Be/NS HMXRB class, 
whose optical counterpart was identified with the star X Persei 
dBraes & Milevll972t(van den Berghll972llBrucato & Kristiaii 
Il972t IWeisskopf et al.lll984b . In this system the companion 
pulsar accretes material from the massive, low-velocity equa- 
torial wind of the Be primary star and emits X-rays as the grav- 
itational energy of the infalling material is converted to X-ray 
luminosity. 

From the optical point-of-view, X Persei (HD 24534, HR 
1209) is a bright and highly variable star (on time scales 
from minutes to years) in optical and infrared (IR) bright- 
ness (V = 6.1-6.8, K = 5.2-6.7) and in emission-line strength 
dMooket alJll974t iRoche et all 1 1997b . When the star is par- 
ticularly bright, its spectrum displays strong emission in H Q 
and other Balmer lines which marks this clearly as a Be- 
type system; when the star is faint, the emission lines dis- 
appear, and it appears to be a norm al early-type star (see: 
Fabregat et al 1992; locheetal 1 19971 and references therein). 



In the past, sometimes it showed also extreme, extended low 
states dNorton et al.lll99ll) ; during these periods the H Q line 



profile changed from emission to absorption, accompani ed by 



a decrease in the I R flux by over a magnitude ( Reynolds et al 
ll992tlRoche etaDll993h . 



Such a variability is commonly supposed to be caused by 
the formation and dissipation of a disk around t he star, as in 
other Be stars. However, between 1994 and 1995 iLiu&Hand 
(2001) observed that the H a EW continued to increase even 
during the optical and near-infrared fading. They found that, 
although both components (line and continuum) arise from the 
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same circumstellar envelope, their emission peaks at different 
evolution stages of the envelope itself, which correspond to 
different distances from the central star. In this scenario, the 
observed behaviour of the two components is attributed to the 
ejection of matter from the Be star as consecutive expanding 
rings. 

The spectral class of the underlying OB s tar has been es- 
timated to be Q9.5 Hi e by ISlettebaS d 1 982b . More recently 



Lvubimkov et al.l (11997) have used both spectroscopic and pho- 



tometric observations of the low-luminosity diskless phase of 
X Persei to infer the characteristics of the visible component; 
their results suggest a star of B0 Ve type, which is likely 
to have a mass of ~ 13-20 Mq and a radius of 5-10 Rq. 
With a projected equatorial velocity vsini of ~ 1200 km s _1 



(Reynold s et al.lll992t iLvubimkov et al.l 119971) the Be star is 



probably seen under a low or intermediate inclination angle. 

Considering UV, optical and I R data taken in diskless and 
near-diskless states limiting et aL ( 1998b modelled the stellar 
photosphere with T e g- = 31000 K and log(g) = 4. Then they 
derived E(B-V) = 0.39 and, assuming R* = 9 Rq, estimated the 
distance to X Persei as 950 ± 200 pc. With a Galactic latitude 
of -17°, a distance estimate of ~ 1 kpc implies that X Persei 
would be ~ 300 pc from the Galactic plane. This is rather far 
off the plane for an early-type star. However, it is within a few 
degrees of the plane of Gould's Belt and the distance estimates 
for X Persei and the inferred extent of Gould's Belt (~ 600 pc; 
Torra etal . 2000) are marginally compatible. 
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The X-ray source has an unusually low luminosity, since 
up to now it has varied on timescales of years between ~ 2.5 
and ~ 15xl0 34 ergs s" 1 ; moreover, while most Be/NS binary 
systems are transient in nature, 4U 0352+309 is a source of 
persistent X-ray emission. This is perhaps due to the relatively 
wide and nearly circular orbit of the system, whose c o mplet e 
ephemeris were determined by iDelgado-Martf et al.1 (120011) . 
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They used measurements made with the All Sky Monitor (ASM) 
on RossiXTE to perform a pulse timing analyses of the 837 s 
pulsations and found strong evidence for the presence of or- 
bital Doppler delays. In this way they inferred that the orbit is 
characterized by a period P or b = 250.3 days, a projected semi- 
major axis of the neutron star a^sim = 454 lt-s, a mass function 
f(M) — I.6IM0 and an eccentricity e = 0.11. The measured 
orbital parameters, together with the known properties of the 
classical Be star X Persei imply a semimajor axis a = 1 .8-2.2 
AU and an orbital inclination i ~ 26-33°. 



massive component and the orbital motion of the neutron star, 
whose trajector y is inclined to t he plane of decretion disk of the 



Pfahl et al. (2002) proposed this object as the prototype for 



a new class of HMXRBs characterised by long orbital peri- 
ods (P or b > 30 days) and low eccentricities. These systems are 
so much wide that no tidal circularisation can have occurred 
after the supernova explosion that created the neutron star, 
therefore all induced eccentricity must be primordial. While 
high velocity kicks are n ecessary to explain the eccentric orbit s 
of most BeX transients (Ivan den Heuvel & van Paradiislll997l) . 
the low-e orbit of X Persei and of the other systems of its class 
must have formed in a supernova explosion without kick. 

Its 8 37 s pulsation period was discovered with the Uhuru 
satellite (Whi te et al.lll976l fl977t iMargon et al.lll977l) and is 
still one of the longest periods of any kno wn accreting pul- 



sar (BildstenetaL|Ll997; Ha berlet al.l ll998). The pulse profile 



of this source is almost sinusoidal. Some authors ( Rob ba et al 



Il996t ICoburn et al.ll2001l) have reported the presence of struc- 
tures in the hardness ratio (4-11 keV)/(1^4 keV), i.e., a sharp 
hardening episode corresponding to the pulse minimum and 
two softening episodes just before and after this hard spike. 

Like many accreting pulsars, the history of the pulse pe- 
riod of 4U 0352+309 shows no evidence for a regular spin- 
up. For the first years after the discovery of pulsations, the 
pulsar exhibited apparently erratic pulse frequency variations 
on a time-scale of a few days, which were superposed on a 
long term trend in which it was spinning up at the rate of 
P/P = —1.5 x 10~ 4 yr _1 . This was followed by a large 
apparent torque reversal around 1978; since then, although 
the spin period measurements have been infrequent, the star 
has been s pinning down at an average rate of P/P = 1.3 x 



1Q- 4 vr" 1 dWhite et all 1 9761: iRobbaet al.l 1 996t Idi Salvo et al 



1998; Delgad o-Martf et al.l2001t and references therein). In all, 
there have been only about two dozen determinations of the 
pulse period of 4U 0352+309, so that its pulse period behavior 
on short timescales (e.g., weeks, months, and even a year) is 
not well documented. 

The transitions between the spin-up and spin-down stages 
occur without any significant variations of the system X-ray 
luminosity and, therefore, cannot be associated with the tran- 
sitions of the neutron star between the accretor and propeller 
states. The long-term spin-down trend argues against the pres- 
ence of a persistent accretion disk in this system; in contrast, 
the assumption of a spherical geometry in the accretion flow 
allows us to interpret the spin period of the neutron star in X 
Persei in terms of the equilibrium period, P^ , provided the 
average (on a time-scale of > 20 yr) value of the wind velocity 
is 350^-00 km s — 1 . The spin-up/spin-down behaviour of the 
source within this scenario can be associated with variations of 
the stellar wind velocity, which are due to the activity of the 



whose trajectory is inclined to the 
Be companion ilkhsanovl2007T) . 



Another unusual aspect of 4U 0352+309 is its X-ray spec- 
trum. Most accreting X-ray pulsar spectra are well described 
by a standard model, in which a simple power law with a pho- 
ton index of ~ 1 is exponentially cut off above ~ 20 keV 
dWhiteet al.l |l983i). The source 4U 0352+309, on the other 



hand, does not follow this standard spectral shape; it is pos- 
sible that this is due to the low mass accretion rate. As a result, 
a variety of different models have been use d to fit the spectrum . 
Using observations taken with RossiXTE ICoburn et al.1 (1200 lb 
found that all these models were improved by the addition of 
a cyclotron resonant scattering feature (CRSF) at ~ 29 keV. In 
particular, the model that best fitted the data is a combination 
of a blackbody plus a power law modified by a CRSF at 28.6 
keV, which implies a magnetic field strength of 2.5(l+z)10 12 G 
in the scattering region (where z is the gravitational redshift). 
In the phase-averaged spectrum there is no evidence of a fluo- 
rescent Fe K-shell emission line in the range 6.4-6.7 keV, and 
the phase-resolved analysis shows that the blackbody and cy- 
clotron line energies are consistent with being constant through 
the pulse. 



Recently. iGrundstrom et al.1 ((2007) have observed that be- 
tween 2000 and 2006 the circumstellar disk grew to near record 
proportions, and concurrently the system X-ray flux dramati- 
cally increased, presumably due to the enhanced mass accre- 
tion from the disk. The H a emission equival ent widths were the 



largest measured over the last few dec ades (IRoche et al.ll 1993 



Piccio ni et al. l200(]Hclarketai]l200lh and they indicated that 
the disk had grown significantly in size. The apparent increase 
in X-ray flux that accompanied the disk expansion confirms 
that the X -ray source is powered b y gas from the Be star disk. 
However, Grundstrom et al. ( 2007b found that the largest H Q 
half-maximum emission radius of the circumstellar disk was 
still much smaller than the separation between the Be star and 
neutron star. Therefore they suggested that tidal forces at peri- 
astron may excite a two-armed spiral in the disk, and that mass 
transfer will occur predominantly after periastron as the gas in 
the arm extension is accreted by the neutron star. The observed 
phase of X-ray maximum, which occurs about one quarter of a 
period after periastron, is consistent with this picture. 

2. Observations and data reduction 

X Persei was observed by XMM-Newton for ~ 31 ks on 
February 2 nd , 2003 . The three EPI C cameras, i.e. o ne pn 



(Striider et al 



2001) and two MOS dTurner et al.ll2001l) . were 



active and operated in Full Frame mode. This gave time res- 
olutions of 73 ms and 2.6 s, respectively. For all of them, the 
Medium thickness filter was used. While in the case of the two 
MOS cameras the effective source exposure time is ~ 3 1 ks, it 
is much lower (^16 ks) for the pn camera, since it was affected 
by the 'counting mode' effect during the observation. 

We used version 7.0 of the XMM-Newton Science Analysis 
System (SAS) to process the event files. After the standard 
pipeline processing, we looked for possible periods of high in- 
strumental background, due to flares of soft protons with ener- 
gies less than a few hundred keV. We found that the first ~ 20 
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ks of the observation were affected by a high soft-proton con- 
tamination. However, since X Persei was detected with a very 
high count-rate (larger than 3 cts s _1 in each camera), much 
higher than that of the background, the soft-protons during the 
bad time intervals have a negligible effect on the source spec- 
tral and timing analysis (their count rate in the source extraction 
area is less than 0.1 cts s _1 ). Hence for our analysis we used 
the data of the whole observation. 

3. Timing analysis 

Due to the very high count rate, the data collected by the three 
EPIC cameras were highly affected by photon pile-up. For this 
reason, we had to exclude the central area of the source PSF 
and to consider only data in an annulus area. 

For the MOS1 camera we selected an extraction corona 
with internal and external radius of 45 and 70", respectively; 
for the MOS2 camera the corresponding parameters are 40 and 
70". In the case of the pn camera the source was imaged close 
to a CCD gap, therefore we were forced to select a very nar- 
row corona, with internal and external radius of 30 and 35", 
respectively. We checked with the SAS task epatplot that, with 
these selections, no event pile-up affected our data, therefore 
we used them to accumulate the source light-curves. The back- 
ground curves were accumulated on large circular areas with no 
sources and radius of 60", 40" and 45" for the MOS1 MOS2 
and pn camera, respectively. 

In Fig. 03 we report the total background-subtracted light 
curves in the 0.15-10 keV energy range. It shows that, dur- 
ing the XMM-Newton observation, there were significant flux 
changes of the source on a few hours time-scale. Variations up 
to ~ 50 % around the average level of ~ 10 cts s -1 are evident. 
Moreover, the source hardness-ratio (HR) between the energy 
bands above and below 2 keV varies up to ~ 30 % around the 
average value of 1 .45 and, over long time-scales, increases with 
the total count rate (Fig. [2j». 



Count Rate (0.15 - 10 keV, c s ') 

Fig. 2. Hardness-ratio variation of X Persei as a function of the 
0.15-10 keV count rate, with a time bin of 1 ks. 

folding analysis using the source events of three cameras. In 
this way we derived a period of 839.3 ± 0.3 s. In Fig. [3] we 
show the folded light curves in the three energy intervals 0.15- 
2, 2-10 and 0. 15-10 keV, together with the folded HR between 
the hard and soft light curves. The pulse profile shows a main 
peak at phase ~ 0-0.1, a secondary peak at phase ~ 0.4 and 
a deep minimum at phase ~ 0.5-0.6. The profile is energy de- 
pendent, because the main peak is at phase ~ and ~ 0.1 in 
the high and low energy ranges, respectively. Moreover, below 
2 keV it shows a plateau in the phase range 0.7-0.9, which 
has no counterpart above 2 keV. This energy dependence is 
confirmed by the HR profile, whose maximum is not aligned 
with the count-rate one. The same plot shows also that there 
is not a simple correlation between the hardness and the total 
count rate: we observe the same HR value at completely differ- 
ent count rate levels, but also very different HR values for the 
same count rate. This indicates that the spectral hardness of X 
Persei does not depend in a simple way on its flux level. 
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Background subtracted light curves of X Persei in the 
ranges 0.15-2, 2-10 and 0.15-10 keV, with a time bin 



obtain a measure of the pulse period, we converted the 
times to the solar system barycenter and performed a 
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Fig. 3. Background subtracted light curves of X Persei in the 
energy ranges 0.15-2, 2-10 and 0.15-10 keV, folded at the 
best-fit period. 

In order to go deeper in our investigation, we have con- 
sidered also a finer energy division by dividing the hard en- 
ergy range in the two sub-ranges 2-4 and 4-10 keV. In Fig. |4] 
we compare the folded soft and hard hardness ratios, i.e. the 
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count-rates ratios in the energy bands (2-4)/(0.15-2) keV and 
(4— 10)/(2—4) keV, respectively. In this way we can show that 
the minimum of the first ratio (and of the total light-curve) co- 
incides with the maximum of the second one. This behaviour 



was already observed by Ro bba et al.l (1 1996) with the data col- 
lected in 1990 by the GINGA satellite. 
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Fig. 4. Hardness ratios of the three light curves of X Persei in 
the energy ranges 0.15-2, 2^1 and 4-10 keV, as a function of 
the pulse phase. The vertical lines indicate the phase intervals 
used for the spectral analysis. 

4. Spectral analysis 

For the source and background spectra we considered the same 
extraction regions we used for the light curves. We accumu- 
lated all the events with pattern range 0-4 (i.e. mono- and bi- 
pixel events) and 0-12 (i.e. from 1 to 4 pixel events) for the 
pn and the two MOS cameras, respectively. 

We generated ad hoc response matrices and ancillary files 
using the SAS tasks rmfgen and arfgen. In order to ensure the 
applicability of the \ 2 statistics, all spectra were rebinned with 
a minimum of 30 counts per bin; they were fitted in the energy 
range 0.3-10 keV using XSPEC 1 1.3.2. All spectral uncertain- 
ties and upper-limits are given at 90 % confidence level for one 
interesting parameter. 

After checking that separate fits of the three spectra gave 
consistent results, we fitted the spectra from the three cam- 
eras simultaneously in order to increase the count statistics and 
to reduce the uncertainties. The fit with an absorbed power- 
law yielded N H = (3.56 ± 0.05) x 10 21 cm~ 2 and pho- 
ton index T = 1.17 ± 0.01, but with large residuals and 
X^/d.o.f. = 2.172/2131. The addition of a black-body com- 
ponent improved the fit quality significantly (Fig. [5}: we ob- 
tained N H = (2.86 ± 0.10) x 10 21 cm" 2 , T = 1.35 ± 0.05 
and kT BB = 1-35 ± 0.03 keV, with xj/d.o.f. = 1.052/2129. 
The emission surface of the thermal component has a radius 
R BB = 392± " m (for d 



1 kpc). The unabsorbed flux in 
fx ~ 1.2 x 10~ 9 erg cm" 
s~\ about 39 % of which is due to the black-body component. 
Assuning a source distance of 1 kpc this flux implies an unab- 
sorbed source luminosity Lx ~ 1.4 x 10 35 erg s _1 . 

We also attempted to fit the spectra by replacing the black- 
body component with alternative emission models, such as 



the energy range 0.3-10 keV is fx ~ 1.2 x 10 9 " "'" 2 
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Fig. 5. Top panel: total spectrum of X Persei with the best-fit 
power-law {dashed line) plus black-body (dotted line) model. 
The spectrum of the MOS1, MOS2 and pn camera are shown 
in black, red and green, respectively. Middle panel: residuals 
(in units of a) between data and model in the case of the single 
power-law. Bottom panel: residuals in the case of the power- 
law plus black-body. 

mekal, thermal bremsstrahlung and broken power-law. In all 
these cases the results were worse than those obtained with the 
black-body model, since the reduced chi-squared was higher, 
the residuals were larger and/or the best-fit parameter values 
were unrealistic. 

The spectra reported in Fig.|5]show that the best-fit model 
still leaves an emission feature at 0.5-0.6 keV. We could fit 
it with a gaussian component at E = 0.58 keV, which has an 
equivalent width of 97 eV. On the other hand, we found no 
narrow Iron K Q lines between 6 and 7 keV, with an upper limit 
on its equivalent width of 0. 1 keV at most. The final best-fit 
parameters are reported in Tab. Q] 

Table 1. Best-fit parameters for the total spectrum of X Persei. 
The model parameters are the galactic interstellar absorption 
Aq, the main power-law component PL, the additional black- 
body component BB and the Gaussian line GL. 



Component Parameter 



Value 



A G 
PL 


r 

Flux @ 1 keV 6 


3 42 +o.3U 

J -^ z -0.08 

1 48+ - 02 
9 48+0.64 


BB 


fcT(keV) 


i -^ z -0.02 




R(m) c 


361±3 


GL 


E Hne (keV) 

rd 

i lme 

EQW (eV) 


9 1O+0.22 
A - LJ -0.32 

"'-15 


d.o.f. 




2126 


xl 




1.021 


a 10 21 cm~ 2 






b 10" 2 ph cm" 


2 s" 1 keV" 1 




c for a source distance of 1 kpc 




d 10" 3 ph cm" 


" 2 s" 1 





5. Phase-resolved spectroscopy 

Due to the hardness-ratio variations observed along the pulse- 
period, we performed a phase-resolved spectroscopy in order 



N. La Palombara & S. Mereghetti: XMM— Newton observation of X Persei 



5 



to study in more detail the source behavior. To this aim, we 
analysed the background subtracted spectra in the three phase 
intervals defined in Fig.S i.e. 0.7-1.05 (phase A), 0.05-0.45 
(phase B), and 0.45-0.7 (phase C). 

The first step was to fit all of them with the best-fit power 
law plus blackbody model of the phase averaged spectrum, 
leaving only the relative normalization factors free to vary. In 
Fig. [6] we report the ratios of the three spectra of each instru- 
ment to these renormalized average models. They show sig- 
nificant residuals for the spectra of phases A and C, which 
are harder and softer than the medium one, respectively. They 
clearly demonstrate the spectral variability as a function of the 
pulse phase. 




ftp 



Energy (keV) 

Fig. 6. Ratios of the spectra corresponding to the phase inter- 
vals shown in Fig. [4] to the average model spectrum. The pa- 
rameters are fixed at the best-fit values for the total spectrum, 
except for the normalizations which are at their best-fit values 
for each phase interval. 

We then fitted the three spectra independently. In all cases, 
the absorbed power-law model was not satisfactory, while the 
addition of a black-body component significantly improved the 
fit. Therefore we used this model for all the spectra, leaving all 
the parameters free to vary: the results are reported in Table [2] 
They show that the absorbing column density is almost costant 
along the phase-pulse, as expected if it is due to the interstellar 
extinction. Moreover, both the power-law and the black-body 
components do not change significantly between phases A and 
B, neither in shape nor in intensity. Instead there is a relevant 
variation in phase C, as the power-law component steepens, 
the black-body temperature increases and the total luminosity 
reduces. 

For completeness, we also looked for the presence, in the 
individual spectra, of the emission gaussian component at 0.58 
keV which was found in the phase-averaged spectrum. We 
found that this line is well defined in the spectra of phases B 
and C, with an equivalent width of ~ 105 and ~ 75 eV, respec- 
tively; on the other hand, it is not well constrained in phase A. 

In order to investigate the relative variations of the two 
components with the period phase, we also simultaneously fit- 
ted the three spectra forcing common values for Nh, T and 
kT BB . In this case we obtained N H = 2.941°;^ x 10 21 cm" 2 , 



Table 2. Best-fit spectral parameters for the phase-resolved 
spectroscopy of X Persei in the case of the independent fit of 
the three spectra. Errors are at a 90 % confidence level for a 
single interesting parameter. Nh, kTBB, and Rbb are measured 
in units of 10 21 cm' 2 , keV, and meters (for a source distance 
of 1 kpc), respectively. 



Spectral 




Phase Interval 




Parameter 


A 


B 


C 


N H 

r 

kT B B 


q+u.iy 

o - zo -0.20 

1 on+0.09 
±.ou_ og 

1 3fi+ 005 

^■<JO_0.05 


2 67 +u ' i4 

z -°'-0.14 
-, on+0.06 
-0.06 
i 9 o+0.05 


z - yy: -0.24 

1 61 +0lr 

i A o+0M 

l-48_ og 


Rbb 


425 


433 


264 


/tot 


1.46 


1.49 


0.82 


/PL 


0.83 


0.97 


0.49 




(57 %) 


(65 %) 


(60 %) 


/bb 


0.63 


0.52 


0.33 




(43 %) 


(35 %) 


(40 %) 


xZ/d.o.f. 


1.028/1465 


1.032/1550 


0.990/964 



a Unabsorbed flux in the energy range 0.3-10 keV, in units of 10 

erg cm" 2 s _1 



r PL = l-40± ;^ and kT BB = l-37±°;°2 keV, with x£/d.o.f. 
= 1.105/3989; the corresponding normalization values are re- 
ported in Table [3] In this interpretation the spectral changes as 
a function of the phase are reproduced by the variations in the 
relative contribution of the two components. From this point of 
view, the power-law component provides the same flux during 
phases A and C, while it has a maximum during phase B. On 
the other hand, the contribution of the black-body component 
progressively decreases from phase A to phase C. 



Table 3. Best-fit values for the black-body and power-law 
normalizations, when the three spectra are fitted simultane- 

121 



ously with common values of Nh (2.94lg;Qg x 10 2i cm" 2 ), 
T (lAOtoH) andkT BB (1.37±g;gl keV). Errors at 
confidence level for a single interesting parameter 



Spectral 
Parameter 



K BB 

/tot 

/PL 

/bb 



Phase Interval 
B 



C 



70 ±3 

465±i| 
1.43 
0.65 

(46 %) 
0.78 

(54 %) 



103 ±3 

38811? 
1.50 
0.96 

(64 %) 
0.54 

(36 %) 



67 ±2 

247±1? 
0.84 
0.62 

(74 %) 
0.22 

(26 %) 



ph cm 



a Intensity of the power-law component in units of 10 

s" 1 keV" 1 at 1 keV 

b Radius of the blackbody component (in metres) for a source distance 
of 1 kpc. 

c Unabsorbed flux in the energy range 0.3-10 keV, in units of 10 

erg cm" 2 s" 1 

In order to confirm that the thermal component varies as a 
function of the rotational phase, we should prove that a con- 
stant black-body component is rejected by the data. To this 
aim, we modified the test model by imposing a common value 
of both kT BB and Rbb for the three spectra, while both T 
and Ipl could vary within them. The resulting best-fit has 
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Nh = (3.19±°;°t) 



kT BB = 1.38±g;gl keVand 
R BB = 357^g m, while the power-law parameters are shown in 
Tab. [4] Also in this case the fit quality is very good (x 2 /d.o.f. 
= 1.135/3989), comparable to the one obtained in the previous 
case. Moreover, in this case we find that during phase C most 
(i.e. ~ 60 %) of the total unabsorbed flux is due to the thermal 
component. 



Table 4. Best-fit values for the black-body and power-law nor- 
malizations, when the three spectra are fitted simultaneously 



!+0.08n 
-0.08 ) 



10 2i cm" 2 ), kT BB 



with common values of Nh ((3.19 

(1.38_q 02 v anu i\BB I u_g 

kpc). Errors are at a 90 % confidence level for a single interest 
ing parameter 



02 keV), and R BB (376±l m for a source distance of 1 



Spectral 




Phase Interval 




Parameter 


A 


B 


C 


r 


1 91+0.0^ 

-0.03 


1 44+U't» 

L -^-0.03 


2 16 +tJ ' U7 






112± 3 


64± 3 


/tot 


1.51 


1.51 


0.85 


/PL 


1.00 


1.00 


0.34 




(66 %) 


(66 %) 


(40 %) 


/bb 


0.51 


0.51 


0.51 




(34 %) 


(34 %) 


(60 %) 



a Intensity of the power-law component in units of 10 3 ph cm 
s" 1 keV" 1 at 1 keV 

b Unabsorbed flux in the energy range 0.3-10 keV, in units of 10~ 

erg cm -2 s _1 



The above results are summarized in Fig.|7J where we show 
the pulse-phase dependence of the black-body temperature, 
the power-law photon index and the unabsorbed flux of the 
two components for both fits. Since they have a similar statis- 
tical quality, we can neither confirm nor deny that the thermal 
component is variable. 

6. Discussion 

In Fig. [8] we report the long term evolution of X Persei lu- 
minosity and spin period since the time of the first obser- 
vations. The flux measured by XMM-Newton corresponds to 
~ 10 35 erg s _1 , which is the highest level of luminosity after 
the end of the outburst in 1975. This result is confirmed also 
by the RossiXTEIASM data, which are reported for compari- 
son: they clearly indicates that the XMM-Newton observation 
was performed just at a luminosity peak of the source, which 
was preceded by a steep increase and followed by a sharp de- 
crease. This luminosity increase should imply a corresponding 
increase in the accretion rate and, in the case of matter trans- 
fer to the NS via an accretion disk, a larger momentum trans- 
fer to the neutron star, therefore we would expect to observe 
a pulsar spin-up. This is not confirmed by the measured pulse 
period, which is greatest than in all the previous observations. 
The pulse period evolution shows that the pulsar spin-down has 
proceeded also during the latest years, in spite of the luminos- 
ity increase. This is an independent evidence that the matter is 
transfered from the Be companion to the NS without forming 
an accretion disk. 



1.4 
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co 
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Fig. 7. Pulse-phase dependence of the black-body tempera- 
ture, the power-law photon index and the unabsorbed flux of 
the two components (in units of 10~ 9 erg cm -2 s _1 ), in the 
case of a common temperature and index (top) and of a com- 
mon black-body (bottom) 

The analysis of the pulse profile between 0.15 and 10 keV 



shows that, contrary to previous findings (e .g.: iRobba et al 
1996t Idi Salvo et alJll998t ICoburn et alJboOlh . at all energies 



the pulse shape is far from being sinusoidal. It has a rich struc- 
ture, with a main and a secondary peak, and varies with the en- 
ergy. As a consequence, the peak of the hardness ratio between 
the count rates above and below 2 keV is not aligned with the 
maximum of the total count rate. It is possible that this change 
in the pulse profile is due to the source luminosity level, which 
was higher than in all the previous detections. The analysis of 
the source status around the pulse minimum provided another 
interesting result: in this phase interval also the hardness-ratio 
(2-4)/(0. 15-2) reaches its minimum value, while the hardness- 
ratio (4-10)/(2-4) has a narrow peak . Th is behaviour was al- 
ready observed by Rob ba et alJ ( 1996b and Coburn et al. ( 2001 ) 
on Ginga and RossiXTE data, respectively, at lower luminosity 
levels (i.e. 2-5 x 10 34 erg s _1 ) and indicates a complex spectral 
evolution. 

From the spectral point-of-view, we estimated a hydro- 
gen column density N H = (3.42±q Qg)xl0 21 cm~ 2 , which is 
more than twice the value of (1.51 ± 0.22) x 10 21 cm~ 2 mea- 
sured by Di Salvo et al. (1998) on BeppoSAX data. However, 
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Fig. 8. Luminosity (top) and pulse-period (bottom) history of 
X Persei from 1972 October 29 to 2007 February 4. The lumi- 
nosities of the top panel are based on the unabs orbed flux in 



the 2- 10 keV and on a source distance of 1 kpc (Teltin g et al 
Il998l) . The data with error-bars are based on RossiXTEIASM 



measurements, which are reported for comparison. 

we checked that this difference is due to the different fit model 
they used (i.e. a high-energy cut-off power-law). For refer- 
ence, we note that optical observations give E(B-V) = 0.395 



dTelting et al.lll998l) : assuming A v = 3.1 E(B-V) and the av- 
erage relation A v = N H X 5.59 x 1 0~ 22 cm" 2 between opti - 



cal extinction and X-ray absorption dPredehl & Schmitll 19951) 
this would predict N H = 2.19 x 10 21 cm -2 . 

It is int eresting to compare our spectral results with those 
obtained bv lCoburn et al.l (1200 lb on RossiXTE data, since they 
fit the source spectrum with a two component model equal to 
our model. They estimated a steeper power-law component 
(T = 1.83 ± 0.03 instead of 1.48 ± 0.02), while the black- 
body component had a comparable temperature (kT = 1 .45 ± 
0.02 keV instead of 1.42± ;^ keV). However, in our case the 
black-body radius is much larger (361 m instead of 130 m): 
this is not surprising, since during the XMM-Newton observa- 
tion X Persei was about 5 times more luminous than during the 
RossiXTE observation. In fact, in spite of this luminosity dif- 



ference, also in the case of the RossiXTE observation ~ 60 % 
of the total flux was due to the power-law component. 

As in previous observations of X Persei, we found no evi- 
dence of an Fe-K line in the source spectrum. We estimated a 
90 % confidence upper limit on the equivalent width of almost 
0.1 keV, i.e. an order of magnitude larger than in other works. 
On the other hand, we found a clue of an emission line at 0.58 
keV, which was never found before. It is rather weak (EQW = 
97 eV) and was not detected in all the pulse-phases, therefore 
its analysis requires further investigations. 

Also Coburn et al. d200ll) performed a phase-resolved 
spectroscopy, even if with a different selection of the phase- 
intervals. Contrary to our findings, in their case both the power- 
law photon-index and the black-body temperature are almost 
constant along the pulse profile, and only the component nor- 
malizations change significantly. However, these results were 
obtained with a much lower count statistics and at a luminosity 
level about 5 times lower than in our case; furthermore, a fixed, 
common value of the hydrogen column density was forced to 
the different spectral phases and only data above 3 keV were 
considered. 

The results we have obtained for X Persei are in agree- 
ment with what has been observed in other two Be/NS bina- 
ries with long puls e period, i.e. in the low luminosi ty system 
RX J0146.9+6121 (lLa Palombara & Mereghettil2006l) and dur- 
ing the low-l uminous quiescent state of the transient sourc e 
3A 0535+262 dOrlandini et al.ll2004UMukheriee & Paulll2005h . 
Also in these cases a thermal excess was detected above the 
main power-law: it has been modeled as a black-body with 
kTeB > 1 keV and Rbb < 1 km, which contributes ~ 30 
% of the total source luminosity. From this point of view, this 
class of low luminosity, long period Be/NS binaries is sig- 
nificantly different from the standard high-luminosity X-ray 
pulsars, whose thermal component is usually characterized by 
temperatures of ~ 0. 1 keV and rad ius of a few hundred km (see 
lLa Palombara & Mereghettill2006l for a review). 



Hickox et al.1 (|2004) showed that, in the low luminousity 



binary pulsars (Lx < 10 36 erg s _1 ), the thermal component 
can be due to either emission by photo-ionized or collisionally 
heated diffuse gas or thermal emission from the surface of the 
neutron star. Based on our results, the first option can be re- 
jected, since we cannot reconcile it with the observed black- 
body spectrum. We therefore favor the interpretation of the 
soft excess in X Persei as thermal emission from the neutron 
star polar cap. If we assume that the source is in the 'accre- 
tor' status, with matter accretion on the NS surface, the black- 
body emitting radius of ~ 360 m is consistent with the ex- 
pected size of the polar cap. In fact, if Mns = 1.4 M and 
Rns = 10 6 cm, the source luminosity of 1.4 x 10 35 erg s _1 

and, adopt- 
= 8.4 x 10 s 

cm (ICampana et al.lll998l) . In this case, based on the relation 
R co z ~ Rns (Rns/R™) ' 5 faickox et al.ll2004 . we would ob- 
tain R co / ~ 345 m, i.e. a value in good agreement with the 
estimated black-body emitting radius. 

Based on the above result, it makes sense to attribute the 
observed thermal component to the emission from the NS po- 
lar caps, as in the case of RX J0146.9+6121 and 3A 0535+262. 



implies an accretion rate M — 7.5 x 10 g s 
ing Bns = 10 12 G, a mag netospheric radius R r 
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If this description is correct, we would expect to observe some 
variability of the soft component along the pulse phase. From 
this point of view the phase-resolved spectroscopy provides no 
conclusive results, since it proves that both a variable and a con- 
stant thermal component can account for the observed spectral 
variability. 

7. Conclusions 

Analyzing a ~ 3 1 ks long XMM-Newton archival observation 
of the Be/NS binary system X Persei, we found that the source 
was detected at a luminosity level of ~ 1.4 x 10 35 erg s _1 , 
which is the highest value ever reported since the large out- 
burst that occured on 1975. This result was independently con- 
firmed by the RossiXTEIASM data, which has been monitoring 
the source during the latest decade. 

The timing analysis confirmed the pulsar spin-down al- 
ready detected in the previous observations. Moreover, it shows 
that the pulse profile has a complex structure and is energy de- 
pendent. 

From the spectral point-of-view, thanks to the high effec- 
tive area of XMM-Newton over a wide energy band, we could 
perform an accurate analysis also at low energies. In this way, 
we detected a data excess above the main power-law compo- 
nent, which contributes for ~ 40 % to the total source lumi- 
nosity. Based on its black-body spectrum, with an emission 
radius of ~ 360 m, this component could be attributed to a 
polar-cap emission. It is interesting to note that the same hy- 
pothesis has been advanced also to describe the spectra of other 
two low-luminosity and long-period Be/NS pulsars, i.e. RX 
J0146.9+6121 and 3A 0535+262 (during its quiescent state). 

We found no evidence of a narrow Iron K a line between 
6 and 7 keV; on the other hand, our data suggest the presence 
of an emission feature at ~ 0.58 keV, whose description and 
nature require further investigations. 

The phase-resolved spectroscopy shows a large spectral 
variability along the pulse period, which can be modeled with 
both a variable and a constant thermal component. This result 
supports but does not confirm the polar-cap origin of this emis- 
sion, therefore further studies will be necessary in the future in 
order to solve this issue. 
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